LocoGSE, a sequence-based genome size estimator for plants.

Extensive research has focused on exploring the range of genome sizes in eukaryotes, with a particular emphasis on land plants, where significant variability has been observed. Accurate estimation of genome size is essential for various research purposes, but existing sequence-based methods have limitations, particularly for low-coverage datasets. In this study, we introduce LocoGSE, a novel genome size estimator designed specifically for low-coverage datasets generated by genome skimming approaches. LocoGSE relies on mapping the reads on single copy consensus proteins without the need for a reference genome assembly. We calibrated LocoGSE using 430 low-coverage Angiosperm genome skimming datasets and compared its performance against other estimators. Our results demonstrate that LocoGSE accurately predicts monoploid genome size even at very low depth of coverage (<1X) and on highly heterozygous samples. Additionally, LocoGSE provides stable estimates across individuals with varying ploidy levels. LocoGSE fills a gap in sequence-based plant genome size estimation by offering a user-friendly and reliable tool that does not rely on high coverage or reference assemblies. We anticipate that LocoGSE will facilitate plant genome size analysis and contribute to evolutionary and ecological studies in the field. Furthermore, at the cost of an initial calibration, LocoGSE can be used in other lineages.

Unzipped genome assemblies of polyploid root-knot nematodes reveal unusual and clade-specific telomeric repeats.

Using long-read sequencing, we assembled and unzipped the polyploid genomes of Meloidogyne incognita, M. javanica and M. arenaria, three of the most devastating plant-parasitic nematodes. We found the canonical nematode telomeric repeat to be missing in these and other Meloidogyne genomes. In addition, we find no evidence for the enzyme telomerase or for orthologs of C. elegans telomere-associated proteins, suggesting alternative lengthening of telomeres. Instead, analyzing our assembled genomes, we identify species-specific composite repeats enriched mostly at one extremity of contigs. These repeats are G-rich, oriented, and transcribed, similarly to canonical telomeric repeats. We confirm them as telomeric using fluorescent in situ hybridization. These repeats are mostly found at one single end of chromosomes in these species. The discovery of unusual and specific complex telomeric repeats opens a plethora of perspectives and highlights the evolutionary diversity of telomeres despite their central roles in senescence, aging, and chromosome integrity.

Tracking population structure and phenology through time using ancient genomes from waterlogged white oak wood.

Whole genome characterizations of crop plants based on ancient DNA have provided unique keys for a better understanding of the evolutionary origins of modern cultivars, the pace and mode of selection underlying their adaptation to new environments and the production of phenotypes of interest. Although forests are among the most biologically rich ecosystems on earth and represent a fundamental resource for human societies, no ancient genome sequences have been generated for trees. This contrasts with the generation of multiple ancient reference genomes for important crops. Here, we sequenced the first ancient tree genomes using two white oak wood remains from Germany dating to the Last Little Ice Age (15th century CE, 7.3× and 4.0×) and one from France dating to the Bronze Age (1700 BCE, 3.4×). We assessed the underlying species and identified one medieval remains as a hybrid between two common oak species (Quercus robur and Q. petraea) and the other two remains as Q. robur. We found that diversity at the global genome level had not changed over time. However, exploratory analyses suggested that a reduction of diversity took place at different time periods. Finally, we determined the timing of leaf unfolding for ancient trees for the first time. The study extends the application of ancient wood beyond the classical proxies of dendroclimatology, dendrochronology, dendroarchaeology and dendroecology, thereby enhancing resolution of inferences on the responses of forest ecosystems to past environmental changes, epidemics and silvicultural practices.

Origin and evolution of the triploid cultivated banana genome.

Most fresh bananas belong to the Cavendish and Gros Michel subgroups. Here, we report chromosome-scale genome assemblies of Cavendish (1.48 Gb) and Gros Michel (1.33 Gb), defining three subgenomes, Ban, Dh and Ze, with Musa acuminata ssp. banksii, malaccensis and zebrina as their major ancestral contributors, respectively. The insertion of repeat sequences in the Fusarium oxysporum f. sp. cubense (Foc) tropical race 4 RGA2 (resistance gene analog 2) promoter was identified in most diploid and triploid bananas. We found that the receptor-like protein (RLP) locus, including Foc race 1-resistant genes, is absent in the Gros Michel Ze subgenome. We identified two NAP (NAC-like, activated by apetala3/pistillata) transcription factor homologs specifically and highly expressed in fruit that directly bind to the promoters of many fruit ripening genes and may be key regulators of fruit ripening. Our genome data should facilitate the breeding and super-domestication of bananas.

Functional diversity of microbial eukaryotes in a meromictic lake: Coupling between metatranscriptomic and a trait-based approach.

The advent of high-throughput sequencing has led to the discovery of a considerable diversity of microbial eukaryotes in aquatic ecosystems, nevertheless, their function and contribution to the trophic food web functioning remain poorly characterized especially in freshwater ecosystems. Based on metabarcoding data obtained from a meromictic lake ecosystem (Pavin, France), we performed a morpho-physio-phenological traits-based approach to infer functional groups of microbial eukaryotes. Metatranscriptomic data were also analysed to assess the metabolic potential of these groups across the diel cycle, size fraction, sampling depth, and periods. Our analysis highlights a huge microbial eukaryotic diversity in the monimolimnion characterized by numerous saprotrophs expressing transcripts related to sulfur and nitrate metabolism as well as dissolved and particulate organic matter degradation. We also describe strong seasonal variations of microbial eukaryotes in the mixolimnion, especially for parasites and mixoplankton. It appears that the water mixing (occurring during spring and autumn) which benefits photosynthetic host communities also promotes parasitic fungi dissemination and over-expression of genes involved in the zoospore phototaxis and stage transition in the parasitic cycle. Mixoplanktonic haptophytes over-expressing photosynthesis-, endocytosis- and phagosome-linked genes under nutrient limitation also suggest that phagotrophy may provide them an advantage over non-phagotrophic phytoplankton.

Imputed genomes of historical horses provide insights into modern breeding.

Historical genomes can provide important insights into recent genomic changes in horses, especially the development of modern breeds. In this study, we characterized 8.7 million genomic variants from a panel of 430 horses from 73 breeds, including newly sequenced genomes from 20 Clydesdales and 10 Shire horses. We used this modern genomic variation to impute the genomes of four historically important horses, consisting of publicly available genomes from 2 Przewalski's horses, 1 Thoroughbred, and a newly sequenced Clydesdale. Using these historical genomes, we identified modern horses with higher genetic similarity to those in the past and unveiled increased inbreeding in recent times. We genotyped variants associated with appearance and behavior to uncover previously unknown characteristics of these important historical horses. Overall, we provide insights into the history of Thoroughbred and Clydesdale breeds and highlight genomic changes in the endangered Przewalski's horse following a century of captive breeding.

Shared pedigree relationships and transmission of unreduced gametes in cultivated banana.

BACKGROUND AND AIMS: Cultivated bananas resulted from inter(sub)specific hybridizations involving Musa species and subspecies (M. acuminata subspecies, M. schizocarpa, M. balbisiana) and the subsequent selection, centuries ago, of hybrids with parthenocarpic, seedless fruits. Cultivars have low fertility and are vegetatively propagated, forming groups of somaclones. Relatively few of them, mainly triploids, are grown on a large scale and characterization of their parental relationships may be useful for breeding strategies. Here we investigate parental relationships and gamete-type contributions among diploid and polyploid banana cultivars. METHODS: We used SNP genotyping data from whole-genome sequencing of 178 banana individuals, including 111 cultivars, 55 wild bananas and 12 synthetic F1 hybrids. We analysed the proportion of SNP sites in accordance with direct parentage with a global statistic and along chromosomes for selected individuals. KEY RESULTS: We characterized parentage relationships for 7 diploid cultivars, 11 triploid cultivars and 1 tetraploid cultivar. Results showed that both diploid and triploid cultivars could have contributed gametes to other banana cultivars. Diploids may have contributed 1x or 2x gametes and triploids 1x to 3x gametes. The Mchare diploid cultivar group, nowadays only found in East Africa, was found as parent of two diploid and eight triploid cultivars. In five of its identified triploid offspring, corresponding to main export or locally popular dessert bananas, Mchare contributed a 2x gamete with full genome restitution without recombination. Analyses of remaining haplotypes in these Mchare offspring suggested ancestral pedigree relationships between different interspecific banana cultivars. CONCLUSIONS: The current cultivated banana resulted from different pathways of formation, with implication of recombined or un-recombined unreduced gametes produced by diploid or triploid cultivars. Identification of dessert banana's parents and the types of gametes they contributed should support the design of breeding strategies.

Telomere DNA length regulation is influenced by seasonal temperature differences in short-lived but not in long-lived reef-building corals.

Telomeres are environment-sensitive regulators of health and aging. Here,we present telomere DNA length analysis of two reef-building coral genera revealing that the long- and short-term water thermal regime is a key driver of between-colony variation across the Pacific Ocean. Notably, there are differences between the two studied genera. The telomere DNA lengths of the short-lived, more stress-sensitive Pocillopora spp. colonies were largely determined by seasonal temperature variation, whereas those of the long-lived, more stress-resistant Porites spp. colonies were insensitive to seasonal patterns, but rather influenced by past thermal anomalies. These results reveal marked differences in telomere DNA length regulation between two evolutionary distant coral genera exhibiting specific life-history traits. We propose that environmentally regulated mechanisms of telomere maintenance are linked to organismal performances, a matter of paramount importance considering the effects of climate change on health.

Ecology of Endozoicomonadaceae in three coral genera across the Pacific Ocean.

Health and resilience of the coral holobiont depend on diverse bacterial communities often dominated by key marine symbionts of the Endozoicomonadaceae family. The factors controlling their distribution and their functional diversity remain, however, poorly known. Here, we study the ecology of Endozoicomonadaceae at an ocean basin-scale by sampling specimens from three coral genera (Pocillopora, Porites, Millepora) on 99 reefs from 32 islands across the Pacific Ocean. The analysis of 2447 metabarcoding and 270 metagenomic samples reveals that each coral genus harbored a distinct new species of Endozoicomonadaceae. These species are composed of nine lineages that have distinct biogeographic patterns. The most common one, found in Pocillopora, appears to be a globally distributed symbiont with distinct metabolic capabilities, including the synthesis of amino acids and vitamins not produced by the host. The other lineages are structured partly by the host genetic lineage in Pocillopora and mainly by the geographic location in Porites. Millepora is more rarely associated to Endozoicomonadaceae. Our results show that different coral genera exhibit distinct strategies of host-Endozoicomonadaceae associations that are defined at the bacteria lineage level.

Host transcriptomic plasticity and photosymbiotic fidelity underpin Pocillopora acclimatization across thermal regimes in the Pacific Ocean.

Heat waves are causing declines in coral reefs globally. Coral thermal responses depend on multiple, interacting drivers, such as past thermal exposure, endosymbiont community composition, and host genotype. This makes the understanding of their relative roles in adaptive and/or plastic responses crucial for anticipating impacts of future warming. Here, we extracted DNA and RNA from 102 Pocillopora colonies collected from 32 sites on 11 islands across the Pacific Ocean to characterize host-photosymbiont fidelity and to investigate patterns of gene expression across a historical thermal gradient. We report high host-photosymbiont fidelity and show that coral and microalgal gene expression respond to different drivers. Differences in photosymbiotic association had only weak impacts on host gene expression, which was more strongly correlated with the historical thermal environment, whereas, photosymbiont gene expression was largely determined by microalgal lineage. Overall, our results reveal a three-tiered strategy of thermal acclimatization in Pocillopora underpinned by host-photosymbiont specificity, host transcriptomic plasticity, and differential photosymbiotic association under extreme warming.

Integrative omics framework for characterization of coral reef ecosystems from the Tara Pacific expedition.

Coral reef science is a fast-growing field propelled by the need to better understand coral health and resilience to devise strategies to slow reef loss resulting from environmental stresses. Key to coral resilience are the symbiotic interactions established within a complex holobiont, i.e. the multipartite assemblages comprising the coral host organism, endosymbiotic dinoflagellates, bacteria, archaea, fungi, and viruses. Tara Pacific is an ambitious project built upon the experience of previous Tara Oceans expeditions, and leveraging state-of-the-art sequencing technologies and analyses to dissect the biodiversity and biocomplexity of the coral holobiont screened across most archipelagos spread throughout the entire Pacific Ocean. Here we detail the Tara Pacific workflow for multi-omics data generation, from sample handling to nucleotide sequence data generation and deposition. This unique multidimensional framework also includes a large amount of concomitant metadata collected side-by-side that provide new assessments of coral reef biodiversity including micro-biodiversity and shape future investigations of coral reef dynamics and their fate in the Anthropocene.

Pervasive tandem duplications and convergent evolution shape coral genomes.

BACKGROUND: Over the last decade, several coral genomes have been sequenced allowing a better understanding of these symbiotic organisms threatened by climate change. Scleractinian corals are reef builders and are central to coral reef ecosystems, providing habitat to a great diversity of species. RESULTS: In the frame of the Tara Pacific expedition, we assemble two coral genomes, Porites lobata and Pocillopora cf. effusa, with vastly improved contiguity that allows us to study the functional organization of these genomes. We annotate their gene catalog and report a relatively higher gene number than that found in other public coral genome sequences, 43,000 and 32,000 genes, respectively. This finding is explained by a high number of tandemly duplicated genes, accounting for almost a third of the predicted genes. We show that these duplicated genes originate from multiple and distinct duplication events throughout the coral lineage. They contribute to the amplification of gene families, mostly related to the immune system and disease resistance, which we suggest to be functionally linked to coral host resilience. CONCLUSIONS: At large, we show the importance of duplicated genes to inform the biology of reef-building corals and provide novel avenues to understand and screen for differences in stress resilience.

Analysis of the P. lividus sea urchin genome highlights contrasting trends of genomic and regulatory evolution in deuterostomes.

Sea urchins are emblematic models in developmental biology and display several characteristics that set them apart from other deuterostomes. To uncover the genomic cues that may underlie these specificities, we generated a chromosome-scale genome assembly for the sea urchin Paracentrotus lividus and an extensive gene expression and epigenetic profiles of its embryonic development. We found that, unlike vertebrates, sea urchins retained ancestral chromosomal linkages but underwent very fast intrachromosomal gene order mixing. We identified a burst of gene duplication in the echinoid lineage and showed that some of these expanded genes have been recruited in novel structures (water vascular system, Aristotle's lantern, and skeletogenic micromere lineage). Finally, we identified gene-regulatory modules conserved between sea urchins and chordates. Our results suggest that gene-regulatory networks controlling development can be conserved despite extensive gene order rearrangement.

Early dispersal of domestic horses into the Great Plains and northern Rockies.

The horse is central to many Indigenous cultures across the American Southwest and the Great Plains. However, when and how horses were first integrated into Indigenous lifeways remain contentious, with extant models derived largely from colonial records. We conducted an interdisciplinary study of an assemblage of historic archaeological horse remains, integrating genomic, isotopic, radiocarbon, and paleopathological evidence. Archaeological and modern North American horses show strong Iberian genetic affinities, with later influx from British sources, but no Viking proximity. Horses rapidly spread from the south into the northern Rockies and central plains by the first half of the 17th century CE, likely through Indigenous exchange networks. They were deeply integrated into Indigenous societies before the arrival of 18th-century European observers, as reflected in herd management, ceremonial practices, and culture.

Oak stands along an elevation gradient have different molecular strategies for regulating bud phenology.

BACKGROUND: Global warming raises serious concerns about the persistence of species and populations locally adapted to their environment, simply because of the shift it produces in their adaptive landscape. For instance, the phenological cycle of tree species may be strongly affected by higher winter temperatures and late frost in spring. Given the variety of ecosystem services they provide, the question of forest tree adaptation has received increasing attention in the scientific community and catalyzed research efforts in ecology, evolutionary biology and functional genomics to study their adaptive capacity to respond to such perturbations. RESULTS: In the present study, we used an elevation gradient in the Pyrenees Mountains to explore the gene expression network underlying dormancy regulation in natural populations of sessile oak stands sampled along an elevation cline and potentially adapted to different climatic conditions mainly driven by temperature. By performing analyses of gene expression in terminal buds we identified genes displaying significant dormancy, elevation or dormancy-by-elevation interaction effects. Our Results highlighted that low- and high-altitude populations have evolved different molecular strategies for minimizing late frost damage and maximizing the growth period, thereby increasing potentially their respective fitness in these contrasting environmental conditions. More particularly, population from high elevation overexpressed genes involved in the inhibition of cell elongation and delaying flowering time while genes involved in cell division and flowering, enabling buds to flush earlier were identified in population from low elevation. CONCLUSION: Our study made it possible to identify key dormancy-by-elevation responsive genes revealing that the stands analyzed in this study have evolved distinct molecular strategies to adapt their bud phenology in response to temperature.

DNA methylation-based profiling of horse archaeological remains for age-at-death and castration.

Age profiling of archaeological bone assemblages can inform on past animal management practices, but is limited by the fragmentary nature of the fossil record and the lack of universal skeletal markers for age. DNA methylation clocks offer new, albeit challenging, alternatives for estimating the age-at-death of ancient individuals. Here, we take advantage of the availability of a DNA methylation clock based on 31,836 CpG sites and dental age markers in horses to assess age predictions in 84 ancient remains. We evaluate our approach using whole-genome sequencing data and develop a capture assay providing reliable estimates for only a fraction of the cost. We also leverage DNA methylation patterns to assess castration practice in the past. Our work opens for a deeper characterization of past husbandry and ritual practices and holds the potential to reveal age mortality profiles in ancient societies, once extended to human remains.

Interspecific introgression patterns reveal the origins of worldwide cultivated bananas in New Guinea.

Hybridizations between Musa species and subspecies, enabled by their transport via human migration, were proposed to have played an important role in banana domestication. We exploited sequencing data of 226 Musaceae accessions, including wild and cultivated accessions, to characterize the inter(sub)specific hybridization pattern that gave rise to cultivated bananas. We identified 11 genetic pools that contributed to cultivars, including two contributors of unknown origin. Informative alleles for each of these genetic pools were pinpointed and used to obtain genome ancestry mosaics of accessions. Diploid and triploid cultivars had genome mosaics involving three up to possibly seven contributors. The simplest mosaics were found for some diploid cultivars from New Guinea, combining three contributors, i.e., banksii and zebrina representing Musa acuminata subspecies and, more unexpectedly, the New Guinean species Musa schizocarpa. Breakpoints of M. schizocarpa introgressions were found to be conserved between New Guinea cultivars and the other analyzed diploid and triploid cultivars. This suggests that plants bearing these M. schizocarpa introgressions were transported from New Guinea and gave rise to currently cultivated bananas. Many cultivars showed contrasted mosaics with predominant ancestry from their geographical origin across Southeast Asia to New Guinea. This revealed that further diversification occurred in different Southeast Asian regions through hybridization with other Musa (sub)species, including two unknown ancestors that we propose to be M. acuminata ssp. halabanensis and a yet to be characterized M. acuminata subspecies. These results highlighted a dynamic crop formation process that was initiated in New Guinea, with subsequent diversification throughout Southeast Asia.

The banana genome hub: a community database for genomics in the Musaceae.

The Banana Genome Hub provides centralized access for genome assemblies, annotations, and the extensive related omics resources available for bananas and banana relatives. A series of tools and unique interfaces are implemented to harness the potential of genomics in bananas, leveraging the power of comparative analysis, while recognizing the differences between datasets. Besides effective genomic tools like BLAST and the JBrowse genome browser, additional interfaces enable advanced gene search and gene family analyses including multiple alignments and phylogenies. A synteny viewer enables the comparison of genome structures between chromosome-scale assemblies. Interfaces for differential expression analyses, metabolic pathways and GO enrichment were also added. A catalogue of variants spanning the banana diversity is made available for exploration, filtering, and export to a wide variety of software. Furthermore, we implemented new ways to graphically explore gene presence-absence in pangenomes as well as genome ancestry mosaics for cultivated bananas. Besides, to guide the community in future sequencing efforts, we provide recommendations for nomenclature of locus tags and a curated list of public genomic resources (assemblies, resequencing, high density genotyping) and upcoming resources-planned, ongoing or not yet public. The Banana Genome Hub aims at supporting the banana scientific community for basic, translational, and applied research and can be accessed at https://banana-genome-hub.southgreen.fr.

The genomic history and global expansion of domestic donkeys.

Donkeys transformed human history as essential beasts of burden for long-distance movement, especially across semi-arid and upland environments. They remain insufficiently studied despite globally expanding and providing key support to low- to middle-income communities. To elucidate their domestication history, we constructed a comprehensive genome panel of 207 modern and 31 ancient donkeys, as well as 15 wild equids. We found a strong phylogeographic structure in modern donkeys that supports a single domestication in Africa ~5000 BCE, followed by further expansions in this continent and Eurasia and ultimately returning to Africa. We uncover a previously unknown genetic lineage in the Levant ~200 BCE, which contributed increasing ancestry toward Asia. Donkey management involved inbreeding and the production of giant bloodlines at a time when mules were essential to the Roman economy and military.

Gene expression and genetic divergence in oak species highlight adaptive genes to soil water constraints.

Drought and waterlogging impede tree growth and may even lead to tree death. Oaks, an emblematic group of tree species, have evolved a range of adaptations to cope with these constraints. The two most widely distributed European species, pedunculate (PO; Quercus robur L.) and sessile oak (SO; Quercus petraea Matt. Lieb), have overlapping ranges, but their respective distribution are highly constrained by local soil conditions. These contrasting ecological preferences between two closely related and frequently hybridizing species constitute a powerful model to explore the functional bases of the adaptive responses in oak. We exposed oak seedlings to waterlogging and drought, conditions typically encountered by the two species in their respective habitats, and studied changes in gene expression in roots using RNA-seq. We identified genes that change in expression between treatments differentially depending on species. These "species × environment"-responsive genes revealed adaptive molecular strategies involving adventitious and lateral root formation, aerenchyma formation in PO, and osmoregulation and ABA regulation in SO. With this experimental design, we also identified genes with different expression between species independently of water conditions imposed. Surprisingly, this category included genes with functions consistent with a role in intrinsic reproductive barriers. Finally, we compared our findings with those for a genome scan of species divergence and found that the expressional candidate genes included numerous highly differentiated genetic markers between the two species. By combining transcriptomic analysis, gene annotation, pathway analyses, as well as genome scan for genetic differentiation among species, we were able to highlight loci likely involved in adaptation of the two species to their respective ecological niches.